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MODELING  OF  ANISOTROPIC  ELECTROMAGNETIC 
REFLECTION  FROM  SEA  ICE 

Kenneth  M.  Golden  and  Stephen  F.  Ackley 


INTRODUCTION 

It  has  been  established  that  there  exists  anisotropy 
in  the  strength  of  the  return  from  impulse  radar  sound¬ 
ings  in  sea  ice  (Campbell  and  Orange  1974,  Kovacs 
and  Morey  1978).  In  this  report  we  quantitatively 
assess  the  role  of  sea  ice  microstructure  in  determining 
the  nature  of  this  anisotropy. 

Anisotropy  and  sea  ice  macrostructure 

The  amplitude  of  a  radar  signal  reflected  from  the 
ice/water  interface  may  depend  markedly  on  the  azi¬ 
muthal  orientation  of  the  linearly  polarizing  antenna  of 
a  radar  system  with  center  frequency  at  1 00  MHz 
(Campbell  and  Orange  1974).  When  the  antenna  is 
positioned  in  the  horizontal  plane  parallel  to  the  sea 
ice  surface  and  rotated  about  a  vertical  axis,  successive 
maxima  and  minima  of  bottom  signal  strength  are 
spaced  at  approximately  90°  intervals  (Fig.  1 ).  In  the 
most  extreme  cases  of  this  anisotropy,  the  signal  ap¬ 
pears  to  be  completely  extinguished.  This  means,  of 
course,  that  the  signal  returned  is  of  such  small  ampli¬ 
tude  that  it  cannot  be  distinguished  from  background 
noise;  i.e.  the  sensitivity  of  the  receiving  instrument 
has  been  surpassed. 

In  addition  to  anisotropy  of  bottom  reflection, 
Kovacs  and  Morey  (1978)  have  observed  variation  in 
surface  reflection,  depending  on  azimuthal  orientation 
of  the  antenna,  when  measurements  were  done  at  a 
center  frequency  of  625  MHz.  The  anisotropy  associ¬ 
ated  with  the  icc/air  interface  at  this  frequency  is,  how¬ 
ever,  much  smaller  than  that  associated  with  the  ice/ 
water  interface  at  100  MHz.  There  does  not  seem  to 
be  a  well-defined  relationship  between  the  directions 
of  minimum  or  maximum  bottom  and  surface  reflec¬ 
tion.  In  some  cases  the  directions  coincide,  in  others 
the  difference  approaches  90°. 


Sea  Wot«r 


Figure  1.  Rotation  of  antenna  to  achieve  reflective  max¬ 
ima  and  minima. 

Anisotropy  of  reflection  from  sea  ice  has  been 
observed  in  both  first-  and  multiyear  ice  with  a  range 
in  thickness  from  25  cm  to  over  2  m.  Freshwater  ice 
has  never  been  observed  to  exhibit  anisotropy  (Camp¬ 
bell  and  Orange  1974).  Anisotropy  is  observed  more 
often  in  first-year  ice  than  in  multiyear  ice  and  is  more 
pronounced  in  the  younger  ice.  Where  first-year  ice  is 
smooth  and  undeformed,  the  azimuthal  orientations  of 
the  antenna  that  correspond  to  maximum  and  minimum 
bottom  reflections  stay  constant  over  distances  of  many 
kilometers.  In  areas  of  broken-up  or  multiyear  pack  ice 
consisting  of  rotated  and  refrozen  floes,  the  antenna 
orientations  vary  drastically  over  distances  of  several 
meters. 

Anisotropy  and  sea  ice  microstructure 

Associated  with  differences  in  sea  ice  macrostructure 
are  differences  in  microstructure.  For  example,  first- 
year  ice  generally  has  a  higher  brine  content  than  multi¬ 
year  ice,  because  the  ice  rejects  brine  with  age.  We  will 
now  review  sea  ire  microstructure  and  its  relation  to 
the  reflective  anisotropy. 

Sea  ice  is  an  inhomogeneous  anisotropic  material 


consisting  primarily  of  pure  ice  surrounding  pockets 
of  brine.  Tiny  spherical  crystals  of  pure  ice  constitute 
the  first  component  of  newly  forming  sea  ice.  These 
spherical  crystals  grow  rapidly  into  disks  and  then  into 
dendritic  plates.  The  plane  of  a  dendritic  plate  coin¬ 
cides  with  the  basal  plane  of  the  ice  crystal  that  com¬ 
prises  it  and  is  perpendicular  to  the  optic  or  c-axis  of 
the  crystal.  In  calm  water,  the  plates  freeze  together 
to  form  a  smooth  skim  of  ice,  with  c-axis  orientation 
primarily  vertical.  In  turbulent  water,  the  plates  are 
tossed  about  and  broken  up  so  that  when  they  congeal 
into  an  ice  surface,  the  c-axis  orientation  is  quite  ran¬ 
dom.  As  growth  continues  downward,  the  c-axes, 
through  a  process  of  geometric  selection,  become  pre¬ 
dominantly  horizontal.  The  layer  within  the  ice  where 
the  c-axis  orientation  changes  from  vertical  or  random 
to  predominantly  horizontal  is  called  the  transition 
zone.  This  zone  is  about  10  to  15  cm  thick  and  occurs 
at  a  depth  anywhere  from  approximately  2  to  60  cm, 
depending  on  growth  conditions  (Weeks  and  Assur 
1967,  Hobbs  1974,  Anderson  and  Weeks  1958,  Weeks 
and  Gow  1979). 

The  ice/water  interface  where  crystal  growth  takes 
place  is  not  smooth.  When  ice  forms  from  sea  water, 
ions  (salts)  are  rejected,  since  salts  are  extremely  in¬ 
soluble  in  ice  (Hobbs  1974).  The  rejection  of  salts 
causes  an  increased  solute  concentration  adjacent  to 
the  slowly  advancing  solid/liquid  interface.  This  in¬ 
creased  concentration  of  salt  and  an  appropriate  tem¬ 
perature  profile  produce  a  zone  of  constitutionally 
supercooled  liquid  below  the  interface  (Weeks  and 
Gow  1979).  Tiller  (1974)  has  indicated  that  the  type 
of  solid/liquid  interface  associated  with  constitution¬ 
ally  supercooled  liquid  is  dendritic  or  cellular,  rather 
than  planar.  Thus,  the  dendritic  plates  perpendicular 
to  the  c-axes  are  elongated  parallel  to  their  basal 
planes  and  extend  into  the  sea  water  below  (Fig.  2). 

As  further  growth  occurs,  the  concentrated  sea  water 
between  the  dendritic  plates  becomes  trapped  and 
makes  up  the  so-called  brine  layers  of  sea  ice.  Since 
the  c-axis  is  perpendicular  to  the  longer  dimensions 
of  the  dendritic  plates,  it  must  also  be  perpendicular 
to  the  longer  dimensions  of  the  brine  layers. 

Immediately  above  the  ice/watcr  interface,  the  ice 
is  relatively  warm  and  contains  a  high  concentration 
of  brine  layers  (Fig.  3a).  t  he  spacing  between  adja¬ 
cent  layers  is  the  thickness  of  a  dendritic  plate.  Far¬ 
ther  up  the  temperature  decreases  and,  consequently, 
the  brine  layers  begin  to  "neck”  (Fig.  3b)  and  change 
into  rows  of  closely  spaced  cylinders  (Fig.  3c),  and 
finally  into  elliptical  cylinders  (Fig.  3d).  Below  the 
transition  zone,  the  layers  or  rows  stand  up  vertically 
with  the  shortest  dimension  horizontal.  Above  the 
transition  zone,  the  arrangement  is  more  random 


Sea  Water 

Figure  2.  Dendritic  plates  at  ice/water  inter¬ 
face,  elongated  parallel  to  their  basal  planes, 
extending  into  the  sea  water  below.  Sea 
water  trapped  between  the  plates  becomes 
the  brine  layers  (from  Weeks  and  Gow  1979). 
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Figure  3.  Brine  layers  (top  view)  near 
the  bottom  of  sea  ice  (a)  begin  to  "neck 
with  decreasing  temperature  farther  up 
in  the  ice  sheet  (b)  and  freeze  out  into 
cylinders  (c)  and  elliptical  cylinders 
(d)  (from  Anderson  and  Weeks  1958). 
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depending  on  growth  conditions.  This  model  of  brine 
structure  is,  of  course,  idealized.  Brine  drainage  influ¬ 
enced  by  gravity  and  temperature  gradient  is  but  one 
factor  that  could  alter  the  idealized  model. 

Below  the  transition  zone  the  c-axes  lie  predom¬ 
inantly  in  the  horizontal  plane.  Accordingly,  layers 
of  brine  in  each  crystal  grain  are  positioned  with  their 
two  longer  dimensions  lying  in  a  vertical  plane  perpen¬ 
dicular  to  the  c-axis,  and  with  their  short  dimension 
oriented  horizontally  and  parallel  to  the  c-axis.  Re¬ 
cent  studies  of  arctic  fast  ice  (Cherepanov  1971 , 

Weeks  and  Gow  1979,  Kovacs  and  Morey  1978)  have 
shown  that  the  c-axes  not  only  have  a  preferred  hori¬ 
zontal  orientation,  but  also  have  a  preferred  azimuthal 
direction  within  the  horizontal  plane.  Cherepanov 
(1971 )  found  the  azimuthal  c-axis  alignment  to  be 
constant  for  hundreds  of  kilometers  in  the  Kara  Sea. 
Weeks  and  Gow  (1979)  and  Kovacs  and  Morey  (1978) 
found  good  positive  correlation  between  the  preferred 
direction  of  c-axis  azimuth  and  “long  term”  current 
direction  beneath  the  growing  ice.  Weeks  and  Gow 
(1979)  proposed  an  explanation  for  the  preferred 
orientation  in  terms  of  these  currents. 

Preferred  azimuthal  orientation  of  c-axes  corres¬ 
ponds  to  preferred  azimuthal  orientation  of  brine  lay¬ 
ers.  Therefore,  in  the  bottom  portions  of  undeformed 
ice,  we  expect  ordered  arrays  of  brine  layers  with  uni¬ 
form  alignment  over  a  large  area.  This  ordered  arrange¬ 
ment  of  c-axes  and  brine  layers  is  apparently  related 
to  the  reflective  anisotropy.  Campbell  and  Orange 
(1974)  and  Kovacs  and  Morey  (1978)  have  found 
that  a  maximum  reflection  from  the  icc/water  inter¬ 
face  occurs  when  the  electromagnetic  wave  is  polar¬ 
ized  so  that  the  direction  of  the  electric  field  parallels 
the  preferred  c-axis  direction,  or  equivalently,  when  it 
is  polarized  so  that  the  electric  field  is  normal  to  the 
principal  brine/ice  interface  (normal  polarization). 

The  minimum  reflection  from  the  ice/water  interface 
occurs  when  the  pulse  is  polarized  so  that  the  direction 
of  the  electric  field  is  perpendicular  to  the  preferred 
c-axis  direction,  or,  equivalently,  the  electric  field  is 
tangential  to  the  principal  brinc/ice  interface  (tan¬ 
gential  polarization).  These  findings  suggest  at  least 
two  sources  for  the  observed  reflective  anisotropy: 

1)  uniform  crystal  arrangement  and  2)  uniform  brine 
layer  arrangement.  Hargreaves  (1977)  has  found  that 
anisotropy  of  reflection  due  to  oriented  crystals  is 
much  smaller  than  that  found  in  sea  ice.  Thus,  we 
reject  uniform  crystal  arrangement  as  the  primary 
cause  of  the  reflective  anisotropy  and  will  offer  a 
theoretical  explanation  for  the  observed  reflective 
anisotropy  of  sea  ice  in  terms  of  the  ordered  arrays  of 
brine  layers. 


A  THEORY  OF  ANISOTROPIC  RADAR 
RETURN  FROM  SEA  ICE 

We  propose  that  the  asymmetrical  geometry  of  the 
brine  layers  causes  an  anisotropy  in  the  penetration  of 
the  impinging  electric  field  into  the  brine  inclusions. 
This  anisotropic  penetration  is  associated  with  an  anis¬ 
otropy  in  the  effective  complex  dielectric  constant  of 
the  sea  ice,  which  determines  the  power  returned  to  the 
receiver. 

Anisotropic  electric  flux  penetration 
into  brine  layers 

To  illustrate  the  proposed  anisotropic  behavior  of 
the  field,  we  take  the  following  values  as  the  dimen¬ 
sions  of  a  representative  bottom  sea  ice  brine  layer. 
Anderson  and  Weeks  (1958)  give  the  horizontal  thick¬ 
ness  b  of  a  typical  brine  layer  as  0.1  mm.  Kovacs  and 
Morey  (1 978)  give  the  vertical  length  c  of  a  brine  layer 
as  5b.  The  horizontal  length  a  of  a  brine  layer  perpen¬ 
dicular  to  the  c-axis,  determined  by  examination  of 
rubbings  from  sea  ice  bottom  cores  (Weeks  and  Assur 
1967),  is  typically  3  mm.  An  idealized  model  of  a 
brine  layer  is  a  non-degenerate  ellipsoid  (a  >  c  >  b)  with 
a  surface  defined  by 

2  2  2 

r+r  +  £_  =  i  (i) 

<z2  7>2  c2 

where  a  -  3  mm,  b  =  0.1  mm,  and  c  =  0.5  mm.  The 
x-axis  is  horizontal  and  perpendicular  to  the  crystal 
optic  axis,  the  y-axis  is  horizontal  and  parallel  to  the 
optic  axis,  and  the  z-axis  is  vertical.  This  analytical 
form  allows  approximation  of  electric  flux  penetration 
into  a  brine  layer. 

Tinga  et  al.  (1973)  derive  an  expression  for  the  elec¬ 
tric  field  E2 k  inside  a  member  of  a  set  of  ellipsoidal  in¬ 
clusions  (brine)  of  complex  dielectric  constant  e2,  rela- 
live  volume  v,  and  uniform  axial  alignment,  dispersed 
in  a  homogeneous  medium  (ice)  of  complex  dielectric 
constant  e, ,  subjected  to  an  initially  uniform  electric 
field  E  directed  along  one  of  the  principal  axes  k  =  a, 
b,  c: 

e2,= _ !l _ 1  (2) 

e,  +rtk(1  -u)(e2-c,) 

where 

„k  - Jl - (3) 

0  (k  +s)\/(a2+s}([,  2+i)(c2+s) 
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is  called  the  depolarization  factor  (Stratton  1941)  of 
the  ellipsoid  along  the  axis  k,  and  depends  only  on  the 
axial  ratios.  Note  that  the  field  inside  the  ellipsoid  is 
uniform.  The  applicability  of  this  equation  to  the  cal¬ 
culation  of  electric  fields  inside  a  brine  layer  (of  larg¬ 
est  dimension  a)  at  100  MHz  depends  primarily  on  the 
condition  that  a  «  X,  which  is  satisfied  for  the  wave¬ 
length  X  =  1 .7  m  in  pure  ice.  In  addition,  since  brine 
volumes  (relative  volume  of  brine  to  ice)  near  the  bot¬ 
tom  of  sea  ice  reach  0.4  and  the  electric  field  inside  a 
brine  layer  is  influenced  by  the  dipole  fields  of  its 
neighbors,  the  derivation  of  the  equation  should  in¬ 
clude  interaction  effects.  In  the  self-consistent  multi¬ 
phase  dielectric  mixture  theory  developed  by  Tinga 
et  al.  (1973),  an  extra  boundary  has  been  introduced 
between  the  inclusions  to  account  for  these  effects  to 
first  order. 

We  assume  the  electric  fields  arc  of  tb-  urn  F(t) 

=  ZT0  e''t‘J,  and  F2 k(f)  =  F 2k0 e',*w,+A*  i  no  spa¬ 
tial  dependence,  where  j  =  ^-1 ,  to  =  2rrx  10a  rad/s 
and  <5  is  the  phase  difference  between  fields.  We 
evaluate  eq  2  for  the  two  cases  of  interest,  k  =  a  and 
k  -  b,  where  na  =  9.67x  I0  J  and  nb  =  8.39x  10"'  are 
calculated  by  numerical  integration.  We  take  e,  =  t’| 
+je'{  -  3.17  +  /0.013  (Johari  and  Charette  1975)  and 
e2  =  e'2+je 2  =  80  +  / 1000  (Vant  1976),  and  set  v  -  t>b 
=  0.29,  where  ub  =  0.29  is  the  brine  volume  at  5  cm 
from  the  bottom  of  a  piece  of  arctic  sea  ice,  as  ob¬ 
tained  from  a  polynomial  fit  of  brine  volume  data. 


normal  component  of  a  harmonically  time-dependent 
electric  field  at  the  interface  between  two  media, 


(6) 


Everywhere  on  the  brine/icc  interface  \ht  normal  com¬ 
ponent  of  the  field  in  the  brine  must  be  very  small 
compared  to  that  in  ice,  |f2n|  =  (3.1 6x  10'3)|£ln|  . 

During  normal  polarization,  the  surface  of  the  brine 
layer  is  primarily  normal  to  the  applied  field  so  that 
the  uniform  field  inside  the  brine  layer  must  be  of  a 
very  small  magnitude.  Physically,  this  exclusion  arises 
from  a  large  depolarizing  field  within  the  brine  layer 
created  by  the  buildup  of  free  (ionic)  charges  and  bound 
polarization  charges  (from  water  dipole  rotations)  on  the 
principal  interface  so  that  the  net  field  inside  the  brine 
layer  is  very  small.  The  existence  of  these  effects  is 
supported  by  the  fact  that  the  relaxation  times  in  brine 
for  ionic  conduction  Tcond=  e2/ coe’^l  .3x  1CT1  0  s 
(Stratton  1941)  and  water  dipole  rotations  Twatcr  *= 

2x  XT' 1  s  (Vant  1976)  arc  both  significantly  shorter  than 
r/2  =  5x  ICT1’  s,  where  r  is  the  temporal  period  of  a  100- 
MHz  wave. 


Case  2.  Tangential  polarization, 
applied  field  F  -  F ,  ||u 

Taking  the  real  part  of  /.'2  to  be  the  physical  scalar 
field,  we  obtain  for  tangential  polarization  from  cq  2, 


Case  I .  Normal  polarization, 
applied  field  F.  -  /  „  ||h 

Taking  the  real  part  of  X2b  =  1 1 ^-2b0  |e',*u’,  +  >l*  to  be 
the  physical  scalar  field  ((ackson  1975  ,  we  obtain 
from  eq  2 

Re(£  ,b)  =  (4.49x10^)  Re(En) 


Re(E2a)  =  (1.90x10',)Re(E1) 

(3.54x  10 ’)  lm(£t)  (7) 

with  h  =  tan ’(3.54x  1 0  ’  / 1  .90x  10  ' )  =61.8°.  The 
moduli  of  the  complex  fields,  or  equivalently  the  am¬ 
plitudes  of  the  physical  fields,  satisfy 


-(5.20x  10  ’)  lm(/:n)  (4) 


where  the  mixing  of  the  real  and  imaginary  parts  of  F. 
indicates  a  phase  shift,  with  fi  =  tan"1  (5.20x  10'1/ 
4.49x  10"4)  =  85.1°.  The  moduli  of  the  complex 
scalar  fields,  or  equivalently  the  amplitudes  of  the 
physical  fields,  satisfy 


“2b  | 


/4,49x  I0^\ 
\  cosX  / 


=  (5.26xl0-’)|En|.  (5) 


m-(^)k,, 

=  (4.01x  10  *  ))£-,) .  (8) 

Therefore,  there  is  a  significant  penetration  of  nearly 
half  the  applied  field  into  the  inclusion  for  tangential 
polarization.  Mathematically,  the  penetration  arises 
from  the  boundary  condition  that  the  tangential  com¬ 
ponent  of  the  electric  field  must  be  continuous  across 
the  brine/icc  interface, 

2«-  <9> 


Therefore,  the  electric  field  is  excluded  from  the  brine 
layer  for  normal  polarization.  Mathematically,  this 
exclusion  arises  from  the  boundary  condition  on  the 
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During  tangential  polarization,  the  surface  of  the  brine 
layer  is  primarily  tangential  to  the  applied  field,  so  that 


a  significant  amount  of  flux  penetrates.  Physically, 
penetration,  or  more  accurately  the  absence  of  a  large 
depolarizing  field,  occurs  because  the  field  is  directed 
along  the  interface  over  most  of  the  surface  causing 
adjacent  polarization  charges  to  cancel  and  free  charges 
to  travel  (though  certainly  not  far)  rather  than  accumu¬ 
late.  The  field  cannot  penetrate  completely,  however, 
since  there  is  a  small  depolarizing  field  created  by  the 
buildup  of  charges  at  the  ends  of  the  brine  layer. 

Implications  of  normal  exclusion,  tangential 
penetration,  and  brine  layer  geometry  for 
dielectric  behavior  of  sea  ice 

Absorption 

Let  S  be  the  surface  of^an  (ellipsoidal)  brine  layer 
V  and  f,  H,  ?2,  H2,  E>2,  B2,  and  J2  be  the  physical 
fields  on  S  and  in  the Jarine  layer.  Then  we  assume 
that  ?  =  ?0  cos(cjz),  E2  =  £20  cos(cj t+S)  and  H2  = 

Wjo  cos(tur+i;)  with  appropriate  dependences  for 
D2,  j2,  and  B2,  where,  in  general,  S  £.  Maxwell's 
equations  may  be  manipulated  (Lorrain  and  Corson 
1970)  to  yield  a  statement  of  energy  conservation  (or 
continuity  of  energy  flow), 


v 


where  n  is  an  outward  unit  normal  vector.  The  rate  at 
which  energy  (represented  by  the  Poynting  vector) 
flows  into  the  brine  layer  is  equal  to  the  increase  in 
electromagnetic  field  energy  stored  in  the  brine  layer 
per  unit  time  plus  the  energy  removed  by  conversion 
to  heat  (conduction  losses)  per  unit  time.  If  one  as¬ 
sumes  that  the  brine  is  linear,  isotropic,  ohmic  and 
non-magnetic,  and  that  the  fields  are  uniform  inside 
the  brine,  then  eq  10  becomes 

~  fixH-n  da  =  {e'2e0E2  +PqH2)  +92^2] 

=  -ojl/|e'2e0£2o  cos(wf+5)  sin(cof +5) 

+  PqH20  cos (wf+£)  sin(wf +£) 

+  <?2^20  cos2(ujf+A)]  (11) 

where  p0  =  4trx  10"7  N/A2  is  the  permeability  of  free 
space,  e0  =  8.854x  10~12  C2/N  m2  is  the  permittivity 
of  free  space,  and  g2  is  the  conductivity  of  brine  in 
mho/m.  If  this  equation  is  divided  by  V  and  time- 


averaged  over  an  integral  number  of  periods,  then  the 
first  two  terms  on  the  right  side  vanish  since  sine  and 
cosine  are  orthogonal,  and  we  are  left  with 

W-\g2E220  (12) 

where  we  have  replaced  the  time-averaged  surface  inte¬ 
gral  with  IV,  the  power  flowing  into  the  volume  V, 
measured  in  W/m3.  For  tangential  polarization,  E2Q 
=  4.01  x  10’1  £0;  for  normal  polarization  £20  =  5.26x 
10‘3  Eq.  Therefore,  the  ratio  of  the  power  losses  is 

-  ’/2g2  (4~01  X  1 0-1  )2  E0  jgggl £0  <  VK)  „  oonn 

W"  '/2?2(S.26x  IO-3)2  Eq  V+/ig]El(V/vb) 

(13) 

where  we  assume  that  9,  is  the  conductivity  of  pure 
ice  and  that  the  amplitude  of  the  field  in  the  ice  is 
equal  to  £0. 

In  other  words,  we  expect  a  very  large  difference  in 
power  absorbed  by  bottom  sea  ice  for  the  two  polar¬ 
izations.  The  variation  is  caused  by  anisotropy  of  elec¬ 
tric  flux  penetration  into  the  brine  layers,  which  in  turn 
is  related  to  anisotropy  in  the  amount  of  brine  layer 
surface  situated  normal  to  the  electric  field. 

Because  the  tangentially  polarized  wave  "sees”  (pen¬ 
etrates)  the  brine  layers  but  the  normally  polarized 
wave  does  not,  significant  conduction  currents  are  set 
up  throughout  the  brine  layer  during  tangential  polar¬ 
ization  but  not  during  normal  polarization.  Note  that 
we  are  primarily  considering  conduction  effects, 
since  100  MHz  is  below  the  resonant  frequency  for 
water  dipole  rotations.  Thus,  for  normal  polarization, 
there  is  a  small  imaginary  mixture  dielectric  constant 
(e”  °c  iVn)  with  a  significant  signal  transmitted  through 
the  medium  and  subsequently  returned  to  the  receiver. 
For  tangential  polarization  there  is  a  large  imaginary 
mixture  dielectric  constant  (e"  <x  ft',)  with  higher  atten¬ 
uation  and  much  less  signal  returned. 

Polarizability 

In  addition  to  anisotropic  absorption  of  the  wave, 
the  asymmetry  of  the  brine  inclusions  causes  aniso¬ 
tropy  in  the  polarizability  of  the  sea  ice,  and  conse¬ 
quently  in  the  real  part  of  the  effective  dielectric  con¬ 
stant.  The  dipole  moment,  defined  (Reitz  and  Milford 
1967)  by 

p=f~rdq  (14) 

for  a  brine  layer  with  surface  5  takes  the  form 
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i  -J 


-  J  ~r  [ot  +  op)da 


where  ot  is  the  tree  (ionic)  surface  charge  density,  op 
is  the  bound  (water  dipole  rotation)  polarization  sur¬ 
face  charge  density,  r  is  a  position  vector,  and  da  is  a 
differential  clement  of  charge.  During  tangential  polar¬ 
ization,  o,  and  Op  contribute  to  the  integral  in  eq  14 
primarily  when  r  is  very  large,  so  that  the  dipole  mo¬ 
ment  p  of  the  brine  layer  is  large.  During  normal  po¬ 
larization,  when  the  applied  field  parallels  the  short 
axis,  op  and  of  are  distributed  more  evenly,  so  that  the 
dipole  moment  of  the  brine  layer  is  smaller.  [See 
van  de  Hulst  (1957)  for  calculations  of  these  dipole 
moments  for  the  two  polarizations.  |  The  dipole  mo¬ 
ments  p  of  the  medium  determine  its  polarization^/5, 
which  is  related  to  the  real  dielectric  constant  by  P  ~ 
£q(c' -  \)T  .  Therefore,  for  tangential  polarization 
there  is  a  larger  real  dielectric  constant  than  for  normal 
polarization. 

We  hope  that  the  preceding  discussion  provides  a 
physical  basis  for  understanding  the  observed  reflective 
anisotropy.  We  feel  that  the  complex  interplay  among 
the  effects  described  and  their  relative  contributions 
to  the  reflective  anisotropy  is  best  examined  with  a 
numerical  model  of  sea  ice  reflections,  which  we  will 
now  present. 


MODELING  OF  ELECTROMAGNETIC 
REFLECTION  FROM  A  STRATIFIED, 
ANISOTROPIC,  INHOMOGENEOUS 
LOSSY  MEDIUM 

We  wish  to  obtain  reflection  profiles  of  sea  ice  of 
thickness  T,  i.c.  graphs  of  power  reflection  coefficients 
(power  received  from  a  given  depth/power  transmitted) 
vs  depth,  for  the  two  polarizations  of  interest.  The  sea 
ice  is  assumed  to  be  composed  of  layers  of  uniform 
thickness  d ,  which  we  shall  take  as  10  cm,  where  the 
number  of  layers  is  C  =  T/d.  (“Layer”  hem  denotes  a 
particular  depth  interval,  c.g.  I  50  to  160  cm.)  To  ob¬ 
tain  a  reflection  profile,  we  calculate  complex  dielec¬ 
tric  constants  c k  ( / )  for  each  layer  /,  where  /'  =  1 , 2,...,tl, 
and  for  each  polarization,  normal  and  tangential. 

These  dielectric  constants  (which  arc  actually  the  first 
two  components  of  the  diagonal  dielectric  tensor)  arc 
used  in  a  calculation  of  interfacial  power  reflection  co¬ 
efficients /?[(/),  where  /  -  1 , 2 . n  and  n  =  (1+ 1 .  Rk(l) 

denotes  reflection  from  the  air/icc  interface  and  P[(n) 
denotes  reflection  from  the  ice/watcr  interface.  The 
interfacial  power  reflection  coefficients  arc  used  in  a 
calculation  of  bulk  power  reflection  coefficients  flk  (/) 


where  /  =  1,2,  since,  for  example,  what  the  re¬ 
ceiver  measures  as  the  bottom  reflection  is  influenced 
by  the  top  reflection.  These  bulk  power  reflection  co¬ 
efficients  are  then  used  to  calculate  attenuated  power 
reflection  coefficients  /?k  (/),  where  /  =  1,2 This 
calculation  takes  losses  (attenuation)  into  account,  as 
determined  by  the  ek(/').  Finally,  the  power  reflection 
coefficients  Rk (/)  are  calculated  from  the  A?k  (/')  by 
including  the  effect  of  beam  spreading.  The  reflective 
behavior  ol  the  radar  pulse  of  small  width  in  frequency 
space  is  approximated  by  analysis  of  the  center  frequency 
(100-MHz)  component  of  a  plane  wave  pulse. 

Calculation  of  mixture  complex 
dielectric  constants 

Tinga  et  al.  (1973)  derive  an  expression  for  the  com¬ 
plex  dielectric  constant  ck  of  a  mixture  consisting  ot 
ellipsoids  (brine)  of  complex  dielectric  constant  e  >, 
relative  volume  v,  and  uniform  axial  alignment,  dispersed 
in  a  homogeneous  medium  (ice)  of  complex  dielectric 
constant  c,  when  the  mixture  is  subjected  to  an  initi¬ 
ally  uniform  electric  field  directed  along  one  ol  the 
principal  axes A  =  a,  b,  c  of  the  ellipsoids: 


ck  =  c  |  + 


uc,  (e2-C|) 
«k(1-u)(c2-c,)  +  c, 


where «k  is  the  depolarization  factor  of  the  ellipsoid 
along  axis  A'.  Since  we  consider  the  sea  ice  to  be  strati¬ 
fied,  cq  16  takes  the  form 


ck(/)  -e,  + 


* "b ( z )  f  t  (f2-C|  ) 

nk(/)(l  -t’t,(/)f(e2-c,)  +  f  | 


=  ck  (/)  +  /' ek(r)  (171 

where  the  brine  volume  i>b(/)  and  depolarization  factors 
/7 k ( / )  of  the  ellipsoids  arc  varied  for  each  layei  r  to 
account  for  the  change  in  physical  properties  of  the  ice 
with  depth,  and  A  denotes  polarization,  either  tangen¬ 
tial  r  or  normal  n.  Equation  1  7  is  evaluated  using  ap¬ 
propriate  frequency-dependent  values  tor  and  c2. 

We  assume  that  the  dielectric  constants  of  the  ice  and 
brine,  Cj  and  e2,  arc  constant  with  depth.  See  Vant 
(1 976)  for  a  discussion  of  the  applicability  of  cq  1 7 
to  time-varying  fields. 

The  brine  volume  uh(/')  is  varied  according  to  poly¬ 
nomial  fits  of  data  calculated  Irom  salinity  and  tempera¬ 
ture  profiles  obtained  from  arctic  fast  ice  by  W.F  .  Weeks 
and  A.J.  Gow.*  Data  are  considered  from  Cape  Krusen- 
stern,  Barrow  (Chukchi  Sea)  and  Harrison  Bay.  The 

*W.I  .  Weeks  and  A.|.  Cow,  CRRLL,  personal  communication,  1979. 
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b.  Calculated  brine  volume  data  with  polynomial  fit  (/.  de¬ 
notes  depth. 


figure  4.  Salinity ,  temperature  and  calculated  brine  volume  data  for  Cape  Krusenstern. 


salinity  and  temperature  profiles  for  the  three  sites  are 
shown  in  Figures  4a,  5a  and  6a,  and  the  brine  volume 
profiles  are  shown  in  Figures  4b,  5b  and  6b.  For  a 
given  10-cm  layer  /',  the  brine  volume  nb(/)  is  calculated 
using  the  equation  ol  f  ranken stein  and  Garner  (given 
in  Weeks  and  Assur  1967), 


i’b(') 


‘I 

io’ 


-49,185 

If 


a  0.532 
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where  5;  is  the  salinity  of  layer  i  expressed  in  parts  per 
thousand,  (I f  is  the  temperature  ol  layer  i  expressed  in 
degrees  Celsius,  and  uh(/)  is  expressed  as  a  dimension¬ 
less  ratio.  Polynomial  tits  were  taken  to  smooth  the 
discontinuities  imposed  by  the  sampling  procedure,  as 
well  as  to  allow  for  variation  of  the  layer  thickness d. 

Recall  that  a  depolarization  factor  tor  a  particular 
axis  depends  only  on  the  axial  ratios  of  the  ellipsoid. 
Thus,  we  vary  the  nk (/)  with  depth  by  changing  the 
relationships  between  the  axes.  In  the  bottom  10  cm, 
we  take  the  relationships  to  be  determined  by  a:h:c  7 
30:1 :5,  so  as  to  model  the  well-defined  brine  layers. 
Throughout  the  Ice  sheet  we  keep  b  constant  but  vary 
a  and  t  linearly  so  that  a:b:c  v  1 :  1 :0.5  at  the  top  of 
the  ice  sheet.  The  linear  change  in  the  ratio a/b  is  to 
simulate  the  necking  and  subsequent  freezing-out  ol 


the  brine  layers  (Fig.  3)  with  a  limit  of  randomly  dis¬ 
persed  tubules  (a  =  b)  at  the  top.  The  vertical  (<  )  dimen¬ 
sion  is  decreased  to  a  value  less  than  a  =  b  at  the  top  to 
indicate  that  the  longest  dimensions  of  the  randomly 
dispersed  tubules  have  a  tendency  to  be  horizontal. 
Therefore,  at  the  top  ol  the  ice  there  arc  fat,  round 
disks  whose  shortest  dimension  is  vertical,  and  at  the 
bottom  there  are  flat,  elongated  disks  whose  shortest 
dimension  is  horizontal.  At  present,  very  little  infor¬ 
mation  is  available  concerning  the  detailed  geometry  ol 
brine  structure.  The  model  we  use  is  a  rough  first  order 
(linear)  approximation  designed  only  to  reflect  the  gross 
features  of  brine  structure:  anisotropy  near  the  bottom 
and  isotropy  near  the  top. 

Calculation  of  interfacial  power 
reflection  coefficients 

We  will  now  calculate  intcrfacial  power  reflection 
coefficients  H[{i)  (power  reflected  by  interlace  //powet 
incident  on  interface  /)  at  each  interface  I  lor  the  two 
polarizations.  The  characteristic  bulk  impedance /k(f) 
of  layer  /  for  polarization  along  axis  k  is  defined  (Ward 
1967)  by 

4(0  09) 

>k(0 
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Depth  (cm) 


Salinity  (%.) 


a.  Salinity  and  temperature. 


Brine  Volume 


b.  Calculated  brine  volume  data  with  polynomial  fit. 

Figure  5.  Salinity,  temperature  and  calculated  brine  vol¬ 
ume  data  for  Barrow  (Chukchi  Sea). 


a.  Salinity  and  temperature. 
Brine  Volume 


b.  Calculated  brine  volume  data  with  polynomial  fit. 

Figure  6.  Salinity,  temperature  and  calculated  brine  vol¬ 
ume  data  for  Harrison  Boy. 
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k 


where  7 k(i)  is  the  propagation  constant  of  layer  /  de¬ 
fined  by 

>k (')  =  Ve Q11QCJ2 e'k(i)  +/p0ojgjd)  •  ^ 

gk(i)  is  the  conductivity  of  layer  /  for  polarization 
along  the  k  axis  of  the  ellipsoids,  and  u  =  2rrx  108 
rad/s.  Since  we  assume  that  the  losses  in  sea  ice  at 
100  MHz  are  primarily  due  to  conduction  effects,  i.e. 
gk(f)  =  e0cuek(/'),  eq  20  becomes 

7k(/')  =  wv/e0d0  ek'('T  •  (21) 


Calculation  of  bulk  power 
reflection  coefficients 

Consider  a  homogeneous  slab  of  pure  ice  (assumed 
lossless)  of  thickness  d  lying  below  a  homogeneous 
half-space  of  air  and  above  a  homogeneous  half-space 
of  sea  water.  Let  a  plane  wave  pulse  of  duration  A t 
traveling  through  air  be  normally  incident  on  interface 
l.  Further  assume  that  A t  «  d/v (,  where  is  the  wave 
velocity  in  ice.  The  primary  pulse  returned  to  the  air 
from  interface  2  has  the  following  bulk  power  reflec¬ 
tion  coefficient  (power  returned  to  air  by  the  inter- 
face/power  incident  on  ice  slab,  without  attenuation): 


The  complex  amplitude  reflection  coefficient  rk(i)  at 
an  interface  /  depends  upon  the  mismatch  of  the  bulk 
impedance  Zk(i-  1)  of  layer  /  -I  above  the  interface 
and  the  bulk  impedance  Zk{i)  of  layer  /  below  the 
interface: 


rk(i)  = 


A('-u+W 


(22) 


On  substitution  of  cq  19  and  21  this  becomes 
r  (/)  _ 


(23) 


«tf(2)=[1-«/(l)|  •  R>(2)  -(1  -R'(\ )|. 

A  secondary  pulse  reflected  back  into  the  air  by  inter¬ 
face  2  has  the  following  bulk  power  reflection  coefficient 

RB( 2)  =  (1-/?'(1)]  •  R\2)  'R\2)  •  [  1  )  | . 

There  are  infinitely  many  of  these  reflections.  For  our 
model  of  a  multi-layered  slab  of  sea  ice  we  neglect  mul¬ 
tiple  reflection  effects  beyond  first  order  and  consider 
only  primary  reflections,  so  that  the  bulk  power  reflec¬ 
tion  coefficient  of  interface  m  for  polarization  k  is 
given  by 


Note  that  this  expression  depends  on  frequency  only 
through  the  dispersive  nature  of  the  dielectric  con¬ 
stants  of  the  constituents  of  the  system,  sea  water,  ice, 
brine,  and  air.  Finally,  the  interfacial  power  reflection 
coefficient  Rk(i)  at  interface  /  is  obtained  from  the 
complex  amplitude  coefficient  through 

RlU)=rkU)*rkU)  (24) 

where  *  denotes  complex  conjugation. 


R%(m)  =  1}  [1-«£(/)|2 R'k(m)  (25) 

i=  1 

for  m  =  2,  3, . . .  ,n,  and 

Rk(m)  =  Rk(m)  (26) 

for  m  -  1 .  Figure  7  gives  a  schematic  of  the  situation. 


Figure  7.  Schematic  diagram  for  bulk  power  reflection  coefficient 
calculation. 
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Calculation  of  attenuated  power 
reflection  coefficients 

The  attenuated  power  reflection  coefficient  R'^(i) 
for  interface  /  (power  returned  to  air  by  interface  // pow¬ 
er  incident  on  ice  slab)  is  calculated  from  R^{i)  by  in¬ 
cluding  losses  from  travel  through  the  sea  ice.  The  am¬ 
plitude  of  a  plane  monochromatic  wave  traveling  one 
way  through  a  lossy  layer  /  of  sea  ice  of  thickness  d 
and  complex  dielectric  constant  e^(/)  is  decreased  by 
the  factor  e  u^,*'/,  where  ct^(i)  is  defined  by  von  Hippcl 
(1954)  as 

°k(/) ;  - (~ '  (v  w  -'|k)  (27) 

where  X0  is  the  free  space  wavelength,  and  k  indicates 
polarization.  One-way  power,  therefore,  decreases  by 
a  factor  [c'ak*'^]2,  which  may  be  termed  a  power 
transmission  coefficient  for  layer  i,  and  wc  define  the 
attenuated  power  reflection  coefficients  R^(m)  by 

*kV)  =  "n  (e“k,'lc'|4/?"(w)  (28) 

i  •  1 

for  m  =  2,3,...  ,n,  and 

K'kV0  =  RlM  (29) 


To  obtain  the  final  power  reflection  coefficients  R^(i), 
we  multiply  each  /?£(/)  by  the  beam  spread  factor  5(/), 

R^)=SU)R^(i)  (33) 


.•ith 


S(i) 


/Jr  (i)  =  /6.83XX2 
Pt  \47t  •  2h  ) 


(34) 


where  h  is  the  distance  Irom  the  antenna  to  interface  /. 
Since  reflection  measurements  are  often  done  with  the 
antenna  resting  on  the  ice,  we  normalize  the  S(r)  so 
that  5(1)  =  1 ;  i.c.  there  is  no  beam  spreading  for  the 
surface  reflection. 

The  above-described  numerical  method  for  obtain¬ 
ing  rel lection  profiles  was  carried  out  on  the  computer 
at  Dartmouth  College,  Hanover,  New  Hampshire.  We 
will  now  present  the  results  of  our  calculations. 


RESULTS 

Anisotropic  bottom  reflections 

Our  calculations  show  anisotropic  bottom  reflections 
from  sea  ice  at  100  MHz.  Figure  8  gives  the  calculated 
reflection  profiles  of  the  two  polarizations  for  the  three 


for  m  -  I . 


Beam  spread 

We  estimate  beam  spreading  effects  with  the  Friis 
transmission  formula  (Kraus  1953),  which  gives  the 
ratio  ol  the  power  in  the  load  of  a  receiving  antenna 
Pr  to  the  power  radiated  by  a  transmitting  antenna 
in  terms  of  the  distance  between  the  antennas r,  the 
wavelength  X,  and  the  directivities  of  the  transmitting 
and  receiving  antenna,  £>,  and  Or,  respectively. 


4rrr 


(30) 


The  di'cctivity  can  be  expressed  as 


t> 

ff(O,0)df  2 


(31) 


where  f(0, <?)  is  the  normalized  power  pattern.  For 
simplicity  we  assume  that  the  transmitting/receiving 
antenna  in  our  study  radiates  energy  uniformlv  over  a 
cone  subtending  a  total  polar  angle  of  90°.  The  direc¬ 
tivity  D  =  Dt  -  Dt  then  becomes 

n  - 1”  .  _ _ *  6.83.  (32) 

/2rd0/"/4  sin DdO 
o  0 


a.  Cape  Krusenstern  (in  the  bottom  layer  vb=  0.1 7). 

Figure  S.  Power  reflection  profile  of  Rfe(i)  for  Cape 
Kusenslern,  Barrow  and  Harrison  Bay.  E,  denotes 
tangential  polarization,  E„  denotes  normal  polariza¬ 
tion. 
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b.  Barrow  (in  the  bottom  layer  ub  =  0. 1 2).  c.  Harrison  Bay  (in  the  bottom  layer  vb  =  0.29). 

Figure  8  (cont  'd). 


Table  1.  Bottom  reflection  anisotropies. 


Site 

0h(  / -C//2) 

K 

Cape  Krusenstern 

0,17 

1  2 

41 

Harrow 

0.12 

1.9 

56 

Harrison  Bay 

0.29 

CO 

800 

sites.  Table  I  gives  the  bottom  layer  brine  volume,  the 
coefficient  of  anisotropy  for  interfacial  bottom  reflec¬ 
tion  K1  =  Rb(n)lR{(n),  and  the  coefficient  of  anisotropy 
tor  net  bottom  reflection  K  =  R n(n)l R t(n)  for  each  site. 
Note  the  correlation  between  bottom  layer  brine  vol¬ 
ume  and  bottom  interfacial  anisotropy  in  Table  I.  This 
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is  not  the  case  for  K  since  net  reflections  are  deter¬ 
mined  by  attenuation  and  reflections  above  the  bottom 
interface. 

Anisotropic  complex  dielectric  constants 

In  our  numerical  model,  reflective  anisotropy  can 
only  exist  with  anisotropic  dielectric  constants.  Figure 
9  shows  the  dielectric  constant  profiles  at  Harrison  Bav 
for  both  normal  and  tangential  polarizations.  The  anis¬ 
otropy  is  most  pronounced  in  the  lower  region  of  the 
ice  where  the  brine  volume  and  alb  ratio  are  high,  i.e. 
where  the  brine  layers  are  numerous  and  well-defined. 
The  wave  is  highly  attenuated  in  these  lower  regions 


o 


(ujo)  qjcteci 


♦ 

v 


b/ 


't 

3 


t 


n 


a.  Normal  polarization.  b.  Tangential  polarization.  Figure  10.  One-way  power  transmission  coeffici¬ 

ents  for  normal  f£n)  and  tangential  (£,)  poiariza - 

Figure  9.  Real  and  imaginary  parts  of  the  complex  mixture  dielectric  constants  at  Harrison  Bay.  tion  at  Harrison  Bay.  Low  transmissions  through 

the  layers  indicate  high  attenuation. 
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a.  a:b:c  =  20:1 :5  at  bottom.  b.  a:b:c  -  50:1:5  at  bottom. 

Figure  1 1.  Power  reflection  coefficient  Rfr(i)  profiles  at  Harrison  Bay. 


because  of  the  very  large  contribution  of  e"  for  tan¬ 
gential  polarization,  as  suggested  previously.  Figure  10 
shows  the  one-way  power  transmission  coefficients 
|e-o(/)j|2  (or  (-Prison  Bay  associated  with  this  atten¬ 
uation.  The  ratio  of  the  imaginary  dielectric  constants 
for  the  bottom  layer  at  this  site,  e"/e'n'  =  47.3/0.027  =*= 
1800  is  roughly  in  agreement  with  the  ratio  of  power 
losses  in  cq  13,  W(/Wn  *  2200,  the  calculation  of 
which  was  based  on  normal  exclusion  and  tangential 
penetration  of  electric  flux. 

We  see  that  the  real  part  of  the  complex  mixture 
dielectric  constant  for  tangential  polarization  at  the 
bottom  of  the  ice  at  Harrison  Bay  has  become  quite 
high  (109),  in  fact  higher  than  that  of  brine  itself  (80). 
The  finite  extent  of  the  brine  layers  causes  them  to 
behave  as  macroscopic  dipoles.  Free  and  polarization 


charges  can  build  up  at  the  ends  of  the  brine  layers, 
rendering  the  sea  ice  highly  polarizable,  more  so  than 
pure  brine  with  no  interfaces. 

The  combination  of  the  large  real  and  imaginary 
parts  of  the  mixture  dielectric  constant  for  tangential 
polarization  cause  the  bulk  impedance  of  the  bottom 
sea  ice  layer  to  be  similar  to  that  of  the  water  below. 
Consequently,  after  downward  propagation,  the  remain¬ 
ing  energy  of  the  tangentially  polarized  wave  tends  to 
be  transmitted,  rather  than  reflected  at  the  ice/watcr 
interface,  as  is  shown  in  Table  1  by  the  K1 .  The  nor¬ 
mally  polarized  wave  encounters  an  impedance  mis¬ 
match  and  is  reflected. 

Sensitivity  of  parameters 

If  one  holds  axes  b  and  c  constant  while  decreasing 
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a,  the  anisotropy  decreases.  Figure  11a  gives  the  re¬ 
flection  profile  for  bottom  a:b:c  =  20: 1 :5.  The  in¬ 
crease  in  bottom  return  for  tangential  polarization  is 
due  to  the  increase  in  surface  area  of  the  ellipsoids 
presented  normally  to  the  field,  which  results  in  less 
electric  flux  penetration  and  less  attenuation,  as  has 
been  previously  shown.  If  the  a  axis  is  increased,  as 
in  Figure  I  lb  where  the  bottom  axial  ratios  are  a:b:c 
~  50: 1 :5,  then  the  return  of  the  tangentially  polarized 
wave  is  decreased  primarily  because  of  higher  attenua¬ 
tion  from  increased  flux  penetration  into  the  brine 
layers.  Note  that  the  bottom  return  for  normal  polar¬ 
ization  is  constant  for  Figures  8c,  I  la  and  1 1b  since 
electric  flux  is  effectively  excluded  from  the  brine 
layers  in  all  cases. 

If  one  holds  axes  a  and  b  constant  while  increasing 
c,  the  anisotropy  decreases;  the  bottom  return  for 
tangential  polarization  increases  while  that  for  normal 
polarization  stays  constant.  The  reasoning  for  this 
effect  is  the  same  as  that  above:  an  increase  in  c  (like 
a  decrease  in  a)  allows  a  larger  depolarizing  field  to  be 
created  in  the  ellipsoid  so  that  less  of  the  wave  will  be 
attenuated.  Decreasing  the  layer  thickness d  smooths 
the  impedance  mismatches  at  interfaces  within  the  sea 
ice  and  thus  decreases  the  magnitudes  of  the  interfacial 
reflection  coefficients. 

Internal  reflection:  the  bumps 

It  is  very  interesting  to  note  the  appearance  of  an 
internal  reflection  located  at  a  depth  of  about  2  773 
(Fig.  8a,  8c,  and  I  la)  for  tangential  polarization.  In 
the  model  this  bump  arises  from  a  superposition  of 
two  effects.  A  large  brine  volume  gradient  and  an  in¬ 
crease  in  axial  ratio  with  depth  first  combine  to  create 
impedance  mismatches  between  layers,  with  resulting 
interfacial  reflections  starting  up  in  the  ice  sheet  and 
extending  downward  to  the  bottom.  When  the  ice 
slab  is  very  lossy  at  the  bottom  (for  example,  when 
a:b:c  =  30 : 1 :5)  the  bottom  reflections  arc  cut  out  by 
attenuation,  leaving  the  bump  higher  in  the  ice  sheet 
(Fig.  8c).  The  position  ot  the  bump  may  be  moved 
lower  in  the  ice  sheet  by  varying  the  axial  ratio  func¬ 
tion  with  depth  so  that  larger  j/b  ratios  are  not  en¬ 
countered  until  the  lowest  portions  ol  the  ice  sheet, 
e.g.  f(z)  -  a  +  d/(z-  T),  where  a  and  0  are  constants 
used  to  fit  the  boundary  conditions  ol  isotropy  near 
the  top  and  maximal  anisotropv  near  the  bottom. 

Kovacs  and  Morey  ( 1 979)  have  postulated  the 
existence  ol  this  bump  in  explaining  how  reflective 
anisotropv  can  coexist  with  travel-time  (and  apparent 
dielectric  constant)  isotropy.  I  hev  hypothesize  that 
the  reflection  associated  with  the  bump  of  tjngcntial 
polarization  (which  they  term  an  upward  shift  of  the 
electromagnetic  boundary)  appears  as  a  small  bottom 


reflection  since  the  higher  dielectric  constant  of  the 
ice  for  tangential  polarization  causes  the  wave  velocity 
to  decrease.  Figure  4  of  Campbell  and  Orange  (1974) 
can  be  interpreted  to  show  return  higher  in  the  ice  sheet 
when  the  antenna  is  oriented  lor  minimum  bottom  re¬ 
flection.  We  leave  this  very  interesting  aspect  ol  the 
reflective  anisotropy  open  to  more  exhaustive  experi¬ 
mentation. 


DISCUSSION 

We  have  assumed  extreme  order  in  the  microstruc¬ 
ture  of  sea  ice.  Effective  axial  ratios  for  brine  inclusions 
at  a  given  depth  or  region  within  sea  ice  are  most  prob¬ 
ably  distributed  about  some  mean,  and  there  is  certainly 
variation  in  the  orientations  of  the  brine  inclusions. 
Flowever,  the  standard  deviations  for  the  preferred  bot¬ 
tom  c-axis  azimuth  given  by  Weeks  and  Gow  (1979) 
arc  between  5°  and  1  5°,  indicating  a  higher  degree  of 
ordering.  Our  purpose  was  to  show  that  if  one  views 
brine  inclusions  as  finite,  anisotropic,  conducting 
bodies  that  arc  small  compared  to  wavelength,  then 
the  observed  reflective  anisotropy  follows  naturally 
from  the  behavior  of  the  electromagnetic  field  around 
and  in  the  bodies,  as  specified  by  Maxwell's  equations. 
To  better  illustrate  mechanisms  determining  reflective 
anisotropy,  we  have  taken  a  simplified  model  of  sea 
ice.  Examination  of  these  limiting  cases  may  render 
the  observed  data  more  understandable. 

For  instance,  in  our  theoretical  explanation  of  re¬ 
flective  anisotropy  we  concern  ourselves  primarily  with 
a  limiting  case  of  brine  structure,  namely  brine  lavers 
(Fig.  3a,  b).  Throughout  much  of  the  ice  sheet,  the 
brine  structure  takes  the  form  of  rows  of  vertical,  cigar¬ 
shaped  inclusions  with  a  top  view  similar  to  that  in 
Figures  3c  and  3d.  However,  even  if  there  is  no  hori¬ 
zontal  anisotropy  in  the  inclusion  (a:b:c  =1:1 :5), 
there  will  be  an  anisotropy  in  electric  flux  penetration, 
as  long  as  the  distance  between  adjacent  inclusions  is 
less  than  the  distance  between  adjacent  rows  ol  inclu¬ 
sions.  During  normal  polarization,  induced  dipoles 
interact  so  that  in  the  region  between  adjacent  inclu¬ 
sions  where  the  electric  field  is  tangential  to  the  surf  ace, 
flux  is  diminished.  Consequently,  penetration  of  the 
field  into  the  "cigars"  for  normal  polarization  is  less 
than  that  for  tangential  polarization,  as  long  as  the 
above  distance  relations  are  kept.  However,  reflective 
anisotropy  associated  with  these  ordered  “cigars”  will 
be  smaller  than  that  associated  with  the  brine  lasers. 

An  initial  impetus  for  the  present  research  was  our 
dissatisfaction  with  the  parallel  plate  waveguide  model 
proposed  bv  Kovacs  and  Morey  (1978)  to  explain  the 
anisotropy  phenomenon.  In  this  model,  the  brine 
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structure  near  the  bottom  of  the  ice  is  assumed  to  be¬ 
have  as  a  parallel  plate  waveguide.  The  frequency  of 
100  MHz  is  well  below  cutoff  for  the  waveguide,  so 
that  attenuation  rather  than  propagation  exists  for 
polarization  parallel  to  the  plates.  The  most  funda¬ 
mental  problem  with  this  theory  is  the  lack  of  evidence 
for  the  existence  of  well-defined  plates  of  the  type  re¬ 
quired  by  their  arguments.  We  feel  that  it  is  more 
natural  to  consider  the  brine  structure  as  consisting 
of  individual,  finite,  anisotropic  bodies  rather  than 
continuous  plates. 

CONCLUSIONS 

The  reflective  anisotropy  observed  in  sea  ice  can  be 
explained  by  consideration  of  the  geometrical  asym¬ 
metries  inherent  in  the  brine  structure  of  sea  ice  with 
a  high  azimuthal  order  of  c-axis  orientation.  Axial 
ratios  of  the  brine  layers  of  the  order  of  30:1 :5  give 
a  reasonable  level  of  power  returned  to  the  ice  surface 
in  the  direction  normal  to  the  brine  layers  (parallel  to 
the  c-axis)  while  significantly  attenuating  the  power 
tangential  to  the  brine  layers  (perpendicular  to  the 
c-axis). 

The  reflective  anisotropy  is  understood  by  consid¬ 
ering  the  amount  of  electric  field  that  penetrates  the 
brine  layers  for  the  tangential  and  normal  polariza¬ 
tions.  Large  electric  field  penetration  in  the  tangential 
case  allows  conduction  effects  to  attenuate  the  wave 
and  reduce  the  power  returned  from  the  icc/water  in¬ 
terface.  Impedance  matching  at  the  ice/water  inter¬ 
face  in  the  tangential  case  also  allows  more  of  the  in¬ 
cident  power  to  penetrate  into  the  sea  water  instead 
of  being  reflected  by  the  interface. 

The  parameters  of  interest  for  the  magnitude  of 
the  reflected  power  in  the  model  calculations  are  the 
brine  volume  profile  with  depth  and  the  axial  ratios 
of  the  brine  inclusions  relative  to  the  polarization  of 
the  wave.  Variation  in  these  parameters  during  model 
calculations  indicated  that  reflection  profiles  would 
be  substantially  changed  by  the  imposed  variations. 
Further  experimental  work  relating  the  observed  re¬ 
flected  power  to  the  details  of  the  sea  ice  microstruc¬ 
ture  (the  axial  ratio  variations  and  distributions  with 
depth,  and  brine  volume  information)  could  be  used 
to  establish  one-to-one  relationships  between  power 
levels  returned  from  certain  depths  and  ice  microstruc¬ 
ture.  If  these  relationships  can  be  established,  radar 
sounding  of  sea  ice  could  then  be  used  as  a  nonde¬ 
structive  "spectroscopic”  tool,  giving  information  on 
the  microstructurc  variations  that  ultimately  control 
the  strength  and  other  important  physical  properties 
of  sev  ice.  At  present,  microstructurc  information 


can  only  be  obtained  by  coring  a  small  sample  of  ice. 

A  methodology  that  could  extend  this  information  to 
line  and  area  measurements  such  as  radar  sounding 
would  be  valuable  in  future  engineering  and  scientific 
investigations  of  sea  ice. 
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